Quantum dynamical screening of the local magnetic moment in Fe-based 

superconductors 
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We have calculated the local magnetic susceptibility of one of the prototypical Fe-based super- 
conductors (LaFeAsO) by means of the local density approximation -I- dynamical mean field theory 
as a function of both (imaginary) time and real frequencies with and without vertex corrections. 
Vertex corrections are essential for obtaining the correct oj-dependence, in particular a pronounced 
low-energy peak at ~ 0.2 eV, which constitutes the hallmark of the dynamical screening of a 
large instantaneous magnetic moment on the Fe atoms. In experiments, however, except for the 
case of x-ray absorption spectroscopy (XAS), the magnetic moment or the susceptibility represent 
typically the average over long time scales. In this respect, the frequency range of typical neutron 
experiments would be too limited to directly estimate the magnitude of the short-time moment. 

PACS numbers: 74.70.Xa, 71.15.Mb, 71.10.Fd, 71.20.Be 



The discovery of high temperature superconductivity 
in iron based pnictidesi not only gave new hope for in- 
creasing the critical temperatures tovifards room temper- 
ature but also posed new challenges for the understand- 
ing of unconventional superconductors, even regarding 
the properties of their normal and magnetic phases. It 
is generally accepted that conventional superconductors 
and also MgB2 are weakly correlated electron systems, 
while cuprate superconductors with a magnetic insu- 
lating ground state, characterized by the opening of a 
Mott-Hubbard gap, are strongly correlated. In con- 
trast, the situation in Fe-based superconductors is much 
less clear. The phase diagram is strikingly similar to 
that of cuprates with an antiferromagnetic parent com- 
pound and a superconducting dome in the doped system. 
However, a major difference from the cuprate physics is 
certainly the metallicity of the antiferromagnetic spin- 
density-wave phase and of the low-doped normal region. 
Also the rather good performance of ab initio (LSDA) 
calculations to describe^ii^ atomic positions and the sym- 
metry of the long-range magnetic order of many Fe-based 
compounds may be interpreted as indicators for the ir- 
relevance of electronic correlations in this compounds. 

In fact, several groups have used weakly correlated 
theories to model iron based superconductors. These 
include the local (spin) density approximation (LDA, 
LSDA), and weak coupling perturbation theory such as 
the random phase approximation (RPA) and the fluctu- 
ation exchange (FLEX) approximation on top of LDA 
bandstructures^^— . The latter show instabilities towards 
magnetism and superconductivity originating from nest- 
ing vectors between different Fermi surface sheets. 

A first hint that a weakly correlated picture is, how- 
ever, insufficient came from angular resolved photoemis- 



sion spectroscopy (ARPES) experimentaii^ showing that 
the LDA bandstructure needs to be rcnormalized by an 
effective mass enhancement of about 2 in LaFeAsO. From 
the theoretical point of view this necessitates the weak 
coupling schemes to start from a correspondingly rcnor- 
malized LDA bandstructure. Such a procedure is partic- 
ular cumbersome for a multi-band system such as iron 
pnictides with five d bands which are all rcnormalized 
somewhat differently and where electronic correlations 
also shift the orbitals relative to each other. In this situ- 
ation a true many-body calculation is almost unavoid- 
able, and indeed several groups have performed such 
calculations. In particular, supplementing the LDA by 
dynamical mean field theory (LDA-fDMFT)^iiS makes 
it possible to describe well the one-particle spectrum 
and renormalization factorsii"— . A second point in- 
dicating that standard band-structure theories do not 
work is the ordered magnetic moment, which is for 
most compounds much larger in LSDA (« 2/^^) than 
in experiment and hardly varies for different Fe-based 
superconductors, in contrast to experiment. In fact, 
the experimental values range from from in 
LaFeAsO (1111 compound)2£i2i via in BaFc2As2 

(122 compound)^ to 2.2^b in FeTe (11 compound)^. 
It is worth also noting here that one generally expects 
the LSDA results to underestimate (rather than over- 
estimate) the magnetic moments for strongly correlated 
electron systems. 

Hence, there is evidently a big moment puzzle with 
theory predicting a large magnetic moments and experi- 
ments measuring smaller ones for most of the Fe-pnictide 
compounds. In a recent Letter—, we put forward a so- 
lution to this puzzle by considering local spin quantum 
fiuctuations, which arc active also above the magnetic 



ordering temperature. In fact, we showed that if one 
considers the local spin-spin correlation function of the 
paramagnetic phase of the Fe-based superconductors on 
very short time scales (fs) a relatively large local moment 
is observed, i.e., comparable to (or even larger than) the 
LSD A prediction of 2/iB for the ordered moment. On 
longer time scales, however, if the electron mobility is 
high enough, this local moment fluctuates very fast so 
that the time-averaged magnetic moment is considerably 
reduced. Such strongly screened moment is the one which 
can become magnetically ordered at low temperatures, 
explaining the reduced size measured in neutron scatter- 
ing experiments. Note that in addition to these temporal 
(but local in space) spin fluctuations, at low-temperature 
also non-local spatial (e.g. antiferromagnetic) spin fluc- 
tuations are expected to become relevant. In this respect, 
only the latter ones are included in the framework of an 
extended Heisenberg model, which has been recently an- 
alyzed in the context of neutron experiments^^ for iron- 
based superconductors. 

From a more theoretical perspective, the origin of the 
physics of large, but dynamically screened, local mag- 
netic moments, as it emerges from our LDA-I-DMFT cal- 
culations, can be related to the presence of several almost 
degenerate moderately correlated Fe-bands at the Fermi 
level. This situation enhances the effects of the Hund's 
coupling interactional, which mainly controls the forma- 
tion of the large local magnetic moment in a still rather 
itinerant metallic background, well fltting to the expres- 
sion "Hund's metals" introduced in this context in Ref. 

m 

Based on our paramagnetic calculations and the itiner- 
ant nature of the antiferromagnetic phase, we suggested 
that the same mechanism not only reduces the local mo- 
ment on long time scales in the paramagnetic phase but 
also the ordered moment in the antiferromagnetic phase, 
corresponding to an average over long time scales. This 
was later confirmed by Yin, Haule and Kotliar-i^iS and, 
independently, by Misawa, Nakamura and ImadaiS^, who 
found a remarkably good agreement with the aforemen- 
tioned experimental moments for different Fe-based su- 
perconductors (see also Ref. 113). The observed mate- 
rial dependence of the ordered magnetic moments comes 
likely from the different degree of itinerancy of the several 
compounds, which results in a different degree of dynami- 
cal screening of the (almost material independent) short- 
time moment. Note that, in the LDA-I-DMFT frame- 
work, another explanation for the magnetic moment mis- 
match has been recently proposed by Hunpyo Lee et al. 
in Ref. H. 

In this paper, we not only extend our earlier studyii 
for different temperatures and provide more details about 
our LDA+DMFT calculation for LaFcAsO, but we also 
present new results for the spectral properties of the (lo- 
cal) spin-correlations functions, aiming to analyze how 
the hallmarks of the dynamically screened local moment 
can be traced in the existing experimental data. 

Specifically, Section U describes how we calculate the 
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FIG. 1. Diagrammatically, the local susceptibility x com- 
puted in DMFT can be expressed in terms of the two (DMFT) 
Green functions, "bubble" contribution (first diagram on the 
r.h.s.), plus vertex (F) corrections in a given particle-hole (in 
our case: spin) channel. 



temporal magnetic moment fiuctuations via the spin-spin 
correlation function. In SectionHH we present new results 
for the temperature dependence of the spin-spin correla- 
tion function. In Section Hill we analytically continue the 
susceptibility to the real frequency axis and discuss the 
importance of vertex corrections. Section IIVI is devoted 
to a detailed comparison and test of the main result of our 
theoretical calculations with experiments. We analyze to 
what extent experimental findings of x-ray absorption 
and neutron scattering spectroscopy are compatible with 
the existence of larger magnetic moments on short (fs) 
time scales, provided an high-frequency extrapolation of 
the experimental data is considered. Finally, we sum- 
marize the main results and our conclusions in Section 
lYl 



I. METHOD AND MODEL 

The most frequently analyzed quantity in 
LDA-fDMFT studies is the spectral func- 
tion, i.e., the single particle Green function 
A{k,uj) = -1/tt Im(G(k,w -f- iO+). This can be 
easily computed in LDA+DMFT and directly compared 
to (inverse-)photoemission experiments. However, the 
DMFT analysis can be also extended, by means of 
an enhanced numerical effort, to the calculation of 
two-particle quantities. For understanding the dynamics 
of the local magnetic moment in the iron-pnictides, 
the quantity of interest is evidently the local spin-spin 
correlation function xioc- In fact, the significant piece 
of information enclosed in xioc has motivated many 
experimental groups to measure or estimate accurately 
Xioc both in the time- and in the frequency- regime: 
As we will discuss more extensively in Sec. IV, these 
experimental estimates represent an important test for 
the relevance of dynamically screened local magnetic 
moments in the physics of iron pnictides, as it emerges 
from our LDA-f DMFT calculations. 

In the paramagnetic case, when the spin-orientation 
can be assumed to be isotropic, the local spin correlation 
function reads 



Xioc.(r) = 2^xIoc.{t) 
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where SI{t) = l/2(n^(T) — nj^(r)) is the z— component of 
the spin operator of the orbital i, expressed in terms of 
the corresponding density operators ^^(t) — c''J{T)d^{T) 
(ft = 1). Note that, in DMFT, we will calculate xioc.(t) 
from the local susceptibility of the converged DMFT im- 
purity model by means of Hirsch-Fye quantum Monte 
Carlo^^. Hence r appearing in Eq. [T] represents here the 
imaginary time. For r > 0, we can write explicitly 

{Si{r)Sm) = l[ ^ «(r)ni(0))-«(r)ni,(0))]. 

(2) 

For the non-interacting case or, similarly, for a given 
quantum Monte-Carlo auxiliary spin configuration^^ {s}, 
this can be expressed in terms of the Green function ma- 

trix G^^{,}(r,0) = (4,(t)c,- ,(0))^^^ as 

{Sl{T)Si{0)) = \[Y1 G%}ir,r)Gi\^^{Q,0)+ 

- (r, 0)G^\,j (0, r) - G:^,y (r, r)G"^ ^^^^ (0, 0)] . 

(3) 

Let us recall that in the present case we work in a local 
basis for which the off-diagonal elements of the Green 
function matrix in the orbital indices vanish within the 
paramagnetic (PM) phase, at least as far as density- 
density interaction are considered. 

After convergence of the DMFT self-consistent loop, 
quantity ^ is measured by means of Monte Carlo sam- 
pling. This procedure represents the most direct way to 
compute local two-particle correlation functions within 
DMFT. In fact, being (numerically) exact at the level 
of the impurity model, it includes automatically all ver- 
tex corrections to the bare-bubble spin susceptibility (see 
Fig. 1, for the diagrammatic illustration of xioc), with- 
out the need of explicitly calculating the local irreducible 
vertex function F. The explicit calculation of F would be 
only needed, in the case one is interested to compute the 
momentum- and frequency dependent spin susceptibility 
x(q, w) in DMFT—. While this provides further possibil- 
ities for the comparison with experiments, for multiband 
systems such as pnictides the calculation of x(q, w) in 
DMFT has been so far possible^ only by introducing 
additional approximations (e.g., neglecting all frequency 
dependencies of the DMFT vertex function F). 

As a last step before turning to the discussion of our 
results, in order to make the comparison with previous 
or new DMFT calculations easier, let us briefly illustrate 
the details of the low-energy Hamiltonian we employed 
for our LDA-hDMFT calculations. 

First of all, the band-structure of LaFeAsO has been 
obtained with the WIEN2k package^, adopting the 
generalized gradient approximation exchange-correlation 
functional introduced by J. P. Perdew^^. As a second 
step, a Wannier projection on a suitable low-energy basis- 
set of maximally localized Wannier functions has been 



performed, using the wien2wannier— and WannierQO^ 
packages. The choice of the most suitable low-energy 
basis-set for DMFT represents clearly an important point 
in the procedure, as it is intrinsically related to the value 
of the Hubbard and Hund's exchange interactions used 
in the calculations, and also to the well-known problem 
of doublc-counting^^i^^ in the LDA-I-DMFT scheme. 

In this respect, let us recall that DMFT calculations 
performed for low-energy model including only the five 
Sd-orbital of Fe run into severe problems connected 
to the different spread of the resulting Wannier func- 
tions. Namely, the difference in the Wannier function 
spreads renders the standard recipes to neglect all 
double counting corrections within the d-manifold 
highly questionable. In order to avoid such problems 
an extension of the basis set for the projection in order 
taking into account the arsenic and oxygen p degrees 
of freedom (resulting in the so called dpp-modcls) is a 
straightforward strategy. 

However, when considering the 1111 systems (in our 
case: LaFeAsO), there is an alternative possibility for 
the construction of a simple low-energy Hamiltonian, re- 
stricted to the d-manifold. The additional approxima- 
tion (first proposed by one of us (R.A.) and H. Ikeda 
m Ref. fe!) can be understood by realizing that in the 
Fc-3(i-manifold, there is specifically one Wannier orbital, 
namely 3z^ — r^, whose localization is much stronger 
w.r.t. to the other four {xy, — y^, xz, yz)^^^!. Hence, 
the local Coulomb interactions for this orbital would be 
different (namely, larger) than for the other ones, making 
no longer justified the neglecting of the double-counting 
correction terms. Fortunately, in the case of LaFeAsO, 
the band mostly connected to the 3z^ — orbital lies 
below the Fermi level and does not contribute to the 
Fermi surface. Hence, one can assume, in fact, that 
only the remaining four bands arc contributing signifi- 
cantly to the response functions under consideration. For 
this reason, we performed the Wannier projection on the 
four remaining Fe3d orbitals only, i.e., xy, — y^, xz, 
yz. This represents the most consistent choice with the 
common assumption of orbital-independent Hubbard and 
Hund's exchange interaction parameters (see Eq. |4] be- 
low), and, consequently, with the neglecting of double 
counting corrections for the d-only low-energy Hamilto- 
nian of LaFeAsO. 

After specifying our four-band low-energy model, we 
consider a density-density type of multi-orbital Hubbard 
interaction, which in second-quantization reads 

4 

-ffint = ^ Un.,^n.,i ^ JSaa')nianjcr'- (4) 

2—1 i^j crcr' 

Note that we set the orbital rotational symmetry relation 
V = U — 2J between the (orbital-independent) Hubbard 
repulsion U for two electrons with opposite spin on the 
same orbital, the off-diagonal repulsion V between dif- 
ferent orbitals i ^ j and the Hund's coupling J. The 
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latter lowers the energy of configurations with parallel 
spins (j = a' . 

The resulting impurity model, associated with DMFT, 
was solved by means of Hirsch-Fye quantum Monte Carlo 
(QMC)^. In agreement with our expectation, the spec- 
tral functions and Fermi surfaces obtained from such a 
four-band model agree much better with the dpp-models 
than the five-band c?-only model does. Moreover, recently 
published data on the spin-spin correlation function ob- 
tained with dpp-modela^ agree well with our previous 
results and, hence, justify the usage of the four band 
model also for the two-particle susceptibilities. 

As mentioned also above, the values of the interac- 
tion parameters depend on the specific choice of the low- 
energy model. Specifically, the Hund's coupling J, stem- 
ming from higher order multipole Coulomb interactions, 
is more robust against the model choice than the Hub- 
bard \J connected to the monopole Coulomb term. Cal- 
culation by means of the constrained random phase ap- 
proximation (cRPA)^ for the five-band model give^ an 
intraorbital Coulomb interaction U of about 2.2 3.3 eV 
and a Hund's rule coupling of about 0.3 0.6 eV. Con- 
sequently, we take an average value of J = 0.45eV while 
we reduce the value of U to 1.8eV in order to take into 
account the screening effects of the eliminated iz^ — 
states and the slightly larger spread of the Wannier func- 
tions (due to the 3z^ — r^-tails). 

Let us stress, at the end of this section that the chosen 
values for the local electronic interaction ([/ = 1.8eV, 
V = U - 23 = 0.9eV, J = 0.45eV) yield the experi- 
mentally observed renormalization of the single particle 
spectra and gives also results in good agreement with the 
dpp-model^. Note that, as described in our previous Let- 
ter, Ref. [13, performing the LDA-I-DMFT calculations 
for an Ising-type Hund's exchange (as in Eq. U]) does not 
lead in the case of the 1111 systems (which are character- 
ized by lower values of the interaction parameters U and 
J) to significant errors, not only for one-particle spec- 
tral functionsii) , but also for the local spin susceptibility 
(inset of Fig. 2 of Ref. [T^). 



II. TEMPERATURE DEPENDENCE 

In the previous Letter we investigated the dichotomy 
between the large instantaneous value of the suscepti- 
bility xiocij) = g^{TrS:,{T)S^{^)) at T = and fixed 
T = 460i^ (/3 = 25 eV~^), which gives us a direct mea- 
sure of the "bare" local moment, and the dynamically 
screened moment at large r, which can be related to the 
moment measured at longer time-scales. We discussed 
the orbital composition of the local spin susceptibility to 
establish that the inter-orbital Hund's coupling is the key 
quantity for the observed behavior. 

In order to provide further insight to the picture of 
magnetic moments driven by local inter-orbital Hund's 
coupling, which get screened at longer time-scales, we 
study here also the temperature dependence of xiocij)- 
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FIG. 2. (color online) Temperature dependence of the 
LDA+DMFT local spin susceptibility as a function of the 
imaginary time r. Inset: Temperature dependence of the lo- 
cal static magnetic susceptibility Xioc('^ = 0), obtained via a 
r-integration. 



In Fig. [2] we show xiocij) for different values of T = 
1//3 ranging from 387K (/? = 30cV-i) to 1160K (/3 ^ 
lOeV""'^). From these results, we note a clear satura- 
tion of the instantaneous spin-spin correlation function 
with decreasing temperature to a value more than twiceas 
large as the corresponding non-interacting one. This 
means that at r = a significant local moment is ac- 
tually formed and it persists at low temperatures. The 
dynamically screened value at t = /3/2, on the other hand, 
decreases evidently upon lowering the temperature. This 
is consistent with what one expects in a metallic system, 
i.e. for a Fermi Liquid. Here, at long enough time scales 
a complete screening takes places at low-T, as a conse- 
quence of the strong fluctuations of the local magnetic 
moment induced by the high electronic mobility (for in- 
stance, in the limiting cases of a non-interacting system 
or of a Fermi Liquid at T = 0, one finds x^J = /?/2) =0). 

Among the various members of the iron-based super- 
conductors, the 1111 family is one of the most itinerant 
ones, therefore the dynamical screening is quite effective. 
As a consequence, the difference between the instanta- 
neous and the moment at r = /3/2 is large and the dy- 
namical screening will be effective at shorter time-scales. 

Moreover, the different temperature dependence of 
XIoc{t = 0) and xioc{j = /3/2) shown in Fig. [2] is also 
refiected in the behavior of the local static susceptibil- 
ity Xioc(w = 0). In fact, our LDA-hDMFT resuUs for 

Xioc{'^ = 0) = /o^ '^'^ Xiocir), reported in the inset 
of Fig. 2, are very weakly dependent on T: We ob- 
serve essentially a strongly renormalized Pauli behav- 
ior. Hence, despite the presence of a large instanta- 
neous magnetic moment, as far as the static moment 
is concerned, we are very far from the Curie behavior 
(Xioc('^ = 0) 7^ const. X ^), which would require, instead, 
that XiociT = 0) and xiociT — l3/2) were almost the same 
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(i.e., a screening time-scale going to oo). A Curie- type 
of behavior is indeed what one obtains for compounds 
with a higher degree of correlation and localization. For 
those, because of the poorer metallic screening, a less 
rapid r-decay of xiociT) is found and the static suscepti- 
bility becomes more Curie-like. 

This is likely the reason why the ordered magnetic 
moment calculated by means of LSDA for some of the 
most correlated ones (e.g., FeTe) is, strangely enough, 
in better agreement with experiments than for the less 
correlated ones. In fact, LSDA is a static theory. This 
means that, even if it cannot reliably describe strongly 
correlated materials, it may work "effectively" better for 
those compounds in which the poor screening flattens the 
XIoc{t), i.e., where it becomes less r-dependent. 

Let us note here, for the sake of completeness, that the 
crossover between Curie-like and (renormalized) Pauli- 
like local spin susceptibility has been nicely illustrated 
by K. Haule, et alM: upon varying the Hund's coupling 
J for a fixed LDA input. For small values of J (close to 
that of the 1111 compounds) xiociT) is quite constant in 
T as in our case, which agrees with the eleven band (dpp) 
calculation for LaFeAsO of Ref . [sgI In contrast, for larger 
values of J the screening of the moments gets less effec- 
tive, and this results in a pronounced 1/T dependence. 
The Hund's coupling J is a crucial parameter in driving 
between these different regimes, in agreement with what 
has been already put forward by W. Yang et al. in the 
context of a x-ray absorption data study^. In a more 
general context, P. Werner, et al^, and, afterwords, L. 
de' Medici, et al^ emphasized that the Hund's coupling 
is a key parameter for tuning correlation effects in multi- 
band systems, with important differences observed de- 
pending on the (integer) filling of the systems^Siil. 



III. FREQUENCY DEPENDENCE 

The calculation of the frequency dependence of the lo- 
cal magnetic spin susceptibility represents an important 
step forward to relate the "Hund's metal" scenario^^^ 
described above with the physics actually observed in 
iron-based superconductors. In fact, in most of the 
experiments (among others neutron experiments, see 
next section) spectral functions (and not time depen- 
dent functions) are measured. At the same time, in 
LDA-I-DMFT the analytic continuation of the local sus- 
ceptibility xioc{t) can be obtained by means of the Max- 
imum Entropy Method^ (MEM) directly from the two- 
particle Green functions of the Anderson impurity model 
associated with the DMFT, namely from the r-dcpcndcnt 
local magnetic spin susceptibility data, as discussed in 
the previous section. In comparison with the analytic 
continuation of the standard one-particle spectral func- 
tions, the numerical effort for accumulating high-enough 
QMC statistic for the MEM is quite larger. Further- 
more, the stability of the analytically continued results 
has been tested by employing different models for the 



MEM, including also a completely featureless one (con- 
stant distribution). 

We start our analysis by considering a rough approx- 
imate expression for Imxioc(w), which is easily obtained 
by neglecting all vertex corrections (F = 0), i.e., by 
considering only the term with the convolution of two 
LDA-hDMFT retarded Green functions G_r(w) (the so- 
called "bubble" term in Fig. [T]). This approximation, 
which is well justified for the evaluation of the optical 
conductivity at the LDA-^DMFT Ievel2ii2,-i^, is expected 
not to work well for the case of spin and charge suscep- 
tibility. We discuss first, however, the "bubble" results, 
because often such easy approximation is taken in the 
theoretical calculations. This preliminary analysis will 
be also useful in the following for evaluating separately 
the effects of the vertex corrections in a situation of prac- 
tical interest, as for the iron-pnictides. 

The spectrum of the local magnetic susceptibility Im 
Xioc(w), computed by considering the "bubble" term 
only, is given by 

2 /> 

Imxioc(c^) = ^ j du:'A{uj')A{uj'^) [/{lu') ~ fico'+iu)] , 

(5) 

where A{uj) = — ilmG/f^w) is the local (k-integrated) 
spectral function obtained by performing the MEM ana- 
lytic continuation of the local LDA-I-DMFT Green func- 
tion, and /(.t) ~ l/(c^^ -t-1) is the Fermi-Dirac distribu- 
tion function. 

The LDA+DMFT resuhs for Imxioc(a;) of LaFcAsO 
at /? = 25eV~"'^ arc shown in Fig. [3] (left panel). The 
sum of the contributions from all four Fe d-Wannier or- 
bitals is shown in the main panel. In the inset, it can be 
seen that the orbital dependence is weak. At the level of 
the "bubble" term, the spectrum of the local magnetic 
susceptibility appears almost featureless, with an unique 
large maximum, located at ~ 2 3eV for all orbitals. 
This has to be related to the fact that (at least) two im- 
portant energy scales, i.e. the LDA bandwidths and the 
local Coulomb interaction in Hint (Eq- HI) are (roughly) 
of this order of magnitude. 

We turn now to evaluate the effects of the inclusion of 
the vertex corrections in the LDA-f-DMFT calculations. 
This corresponds to apply the MEM directly to the two- 
particle local spin susceptibility XIoc(t) as a function of 
imaginary time, without resorting to any expression such 
as Eq. ([S]) in terms of the single-particle Green func- 
tion Gfl(w) only. The results of the (bosonic) MEM for 
XIoc(t) are shown in Fig. [3] (right panel), where the sum 
of all four diagonal contributions (blue line), as well as 
the total sum of all (sixteen) diagonal and off-diagonal 
term (solid thick gray line) are reported. From a com- 
parison of the total intensity of Imx;oc('^), with the bare 
"bubble" data of Fig. [3] (left panel) immediately emerges 
the important role played by the vertex corrections. They 
indeed determine the appearance of a well defined peak 
structure in the spectral functions. Specifically, by look- 
ing at the total intensity of Fig. |3] we note a marked 
low-energy peak, located at w ~ 0.2eV, while a second 
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FIG. 3. Left panel: Frequency dependence of the imaginary part of the local spin susceptibility calculated in LDA+DMFT for 
LaFaAsO at /? = 25eV^^ without vertex corrections ("bare bubble" term, see Eq. ((5|); inset: orbital resolved contributions to 
luixioc{iJj) (note that only the diagonal terms are shown, as the non diagonal are identically zero for the bare "bubble"). Right 
panel: Frequency dependence of the imaginary part of the local spin susceptibility calculated in LDA+DMFT with all vertex 
corrections included: in blue are shown the sum of all four diagonal terms, while the sum of all (sixteen) terms is plotted in 
gray; inset: orbitally resolved plot of the diagonal terms of Im Xioc(i^) after the inclusion of vertex corrections. 



(but much broader) spectral structure is found at larger 
frequencies. 

From the physical point of view, we can look at the 
development of the low-energy peak as the Gngcrprint 
of a local moment formation. Such local moment is, 
however, dynamically screened, as it is witnessed by the 
relatively visible spectral features at higher frequencies, 
and - above all - by the position of the peak itself, at 
ujpeak ~ 150+200 meV. We can better understand this ar- 
gument, by considering a Fourier analysis of Imx(a;), i.e. 
its behavior for real times. The peak of Im%(a;) can be 
directly related to temporal fluctuations of the local spin 
on the Fe atoms with a characteristic frequency ~ uipeak ■ 
These correspond, evidently, to an "effective" decay of 
the local spin moment, when averaging over times larger 
than the typical fluctuating period, whose size can es- 
timated {tfiuct ~ uj^ ^'■^ '-'^ order of ~ 10 fem- 
toseconds (fs). Note that this is larger than the typical 
time-scales (tei) of about 1 fs characterizing the electron 
dynamics in LaFeAsO, i.e. tei ~ (where W ^ AeV is 
the electronic bandwidth, and Z ^ 0.5 the quasiparticle 
renormalization factor), which clearly indicates the rele- 
vance of the local magnetic moment formation even in the 
most metallic family (1111) of the iron-pnictides. At the 
same time, in contrast to the more correlated case of met- 
als in the proximity of a transition to a Mott insulating 
phase (where ujpeak — ^ 0), the local spin fluctuations in 
LaFeAsO occur much faster than the typical time-scales 
of significant experimental probes, such as in the case of 
neutron spectroscopy (see next section for quantitative 
details of the comparison with experiments). 

A second important characteristic sign in Im xi'-^) is 
the importance of the contribution of the off-diagonal 
terms in the spin-correlation matrix, which is reflected 



in the difference between the total and the diagonal part 
of the Imxzoc(w). The significant off-diagonal contribu- 
tion, which stems entirely from vertex correction effects, 
can be easily connected with the predominant role of the 
Hund's exchange interaction in this system. 

Finally, as for the diagonal terms, it is also interesting 
to disentangle our LDA+DMFT data for the local mag- 
netic susceptibility with respect to the four different Fe 
d-orbitals of our model. Contrary to the case of the "bub- 
ble" approximation (inset of the left panel of Fig. [3]), the 
spectra (with vertex corrections included) show visible 
differences among the different orbitals (inset of the right 
panel of Fig. [3]). In particular, one notes that the first 
peak structure at low-energy is much more pronounced 
for the — orbital, i.e., for that whose Wannier func- 
tion lives most in the plane and whose lobes point in the 
direction of the ligand As atoms. One also notes the pres- 
ence of a second peak at energies slightly smaller than the 
second maximum in the total (diagonal) spectrum. On 
the other hand, we note that in the case of the xy orbital 
(which also lives mainly in the plane, but whose lobes 
point in the direction of the neighboring Fe atoms), the 
first peak is considerably weakened. There also appears 
to be a small shift in frequencies which might however be 
an artifact stemming from different peak heights, since 
a combined screening of all orbitals would suggest the 
first peak to be at the same position. The intensity of 
the second maximum is also slightly weakened, and the 
missing spectral weight appears to be shifted towards 
higher frequency, where a third maximum appears (at 
cj 1.7+ 1.8eV). The situation for the degenerate xz, yz 
orbitals is, instead, something in between those of the two 
planar orbitals, with the exception of the position of the 
first peak, which appears moving (very slightly) towards 
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lower energies. The emergence of a stronger low-energy 
peak for the — y^-orbital appears connected to the 
value of the orbital occupations, which are respectively 
nxz = ^yz = 1.06, n^2_y2 = 0.98 and rixy = 0.90. In 
fact, the — orbital is the "closest" to the half-filling 
(ri = 1) condition, and, hence, one expects more evident 
effects of strong-correlations (among which the formation 
of a stronger local magnetic moment). 

The evaluation of the vertex correction effects on the 
LDA-I-DMFT local magnetic susceptibility allows not 
only for a more precise description of the physics charac- 
terizing the different orbitals and of the most-relevant off- 
diagonal contributions, but it also represents a necessary 
step for performing a more quantitative comparison with 
experiments. In fact, the question of how our theoreti- 
cal data really compare with spectroscopic experiments 
in the paramagnetic phase will be addressed extensively 
in the next section. 



IV. COMPARISON WITH EXPERIMENT 

A major outcome of our LDA+DMFT calculation is 
the formation of a large instantaneous (t = 0) local mo- 
ment, mainly driven by the Hund's exchange interac- 
tion J among the four almost degenerate Fe-Sd Wannier 
orbitals at the Fermi level (Sec. |n|. However, as the 
(imaginary) time and frequency dependence (Sec. IIII[) of 
the LDA-I-DMFT local spin susceptibility shows, in the 
case of LaFcAsO such large local moment is efficiently 
screened over large time scales by dynamical fiuctuations. 
When comparing our LDA+DMFT results with experi- 
ments performed in the paramagnetic phase of the Fe- 
based superconductors, hence, the focus should be on (i) 
the magnitude of the local magnetic moment at t = 
and (ii) its dynamical screening over long time scales. 

As for (i), one can extract from our LDA-I-DMFT 
x(t = 0) value quantitative estimates for the instanta- 
neous local magnetic moment rnioc,t=o-, and for the "ef- 
fective" magnitude of the total spin Seff of the Fe sites. 
Specifically, after having extrapolated our LDA-I-DMFT 
data down to T = 200K, and assuming a perfect 
spin-isotropy for the PM phase, we find mioc,t=o = 

9^/{sf) + {S^HW) - V^XiociT = 0) ~ 

which would correspond (m = g^/ Seff {Seff + 1)) to 
an "effective" spin configuration Seff ^ 1.4 for each 
Fe-atom. Note, that Ising and Heisenberg Hund's ex- 
change yield a similar (S*^) in the itinerant regime of 
LaFeAsOii, albeit of course not in more strongly cou- 
pled materials. Such estimates are similar to those 
(jnioc,t=o — 'i-4:fJ.B, Sef f ~ 1.3) obtained in recent 
LDA-I-DMFT calculations^ for an 11-band dpp model 
of LaFeAsO. Even higher values for the local magnetic 
moment were extracted by previous LDA-I-DMFT calcu- 
lations {Seff = 1.8)^, which assume a much larger U 
and J value than that estimated by cRPA for LaFeAsO. 
From the experimental point of view, the detection of 



a local and instantaneous magnetic moment of the Fe- 
atoms would evidently require a very fast experimental 
probe, as, e.g., x-ray absorption spectroscopy (XAS). In 
fact, XAS measurements of Fe ^2,3 edge by T. KroU et 
aJ4^ for LaFeAsO have been fitted by multiplet clus- 
ter calculations, suggesting a high-spin ground state (i.e. 
Seff = 2 even larger than the LDA-I-DMFT predictions). 

It is important, however, to test experimentally also 
the second main aspect of our LDA+DMFT results, i.e., 
the dynamical screening of such large local magnetic mo- 
ments over longer time scales. This is possible by inelastic 
neutron spectroscopy. In fact, neutron experiments have 
been performed for T > T/v, e.g., in the case of CaFe2As2 
(at T = 220K > = 172K)ii and (optimally doped) 
Ba(Fei_:,Co)2As2.-^. 

Inelastic neutron spectroscopy are typically re-casted 
in terms of a q- and oj-dependent dynamical structure 
factor, and this, in turn, can be related to the corre- 
sponding q- and w- dependent magnetic susceptibility. 
Although the crystal field (CF) splitting of iron-pnictidcs 
(e.g., CF ^ 200 — 300 meV) is usually much lower than 
the typical values of 1 — 2eV of other 3d transition metal 
compounds, in the post-processing of the neutron experi- 
mental data for the Fe pnictides the assumption of a com- 
plete quenching of the orbital moment is usually done. 
Hence, the dynamical structure factor is directly related 
to the imaginary part of the (reduced)^ magnetic spin 
susceptibility, i.e., to its spectral function Imx(q, w). 

In our case, we are interested in verifying the presence 
of local magnetic moments. Hence, a Fourier transform 
to real space, i.e., an integral over the momenta q is 
necessary. While neutron scattering data for CaFe2As2 
and Ba(Fei_2:Co)2As2 have been taken only for a given 
set of q values, these data have been used for extract- 
ing the remaining part of the q-dcpendcnce of lmx(q, uj) 
via a fitting proccdur o^^'^^i^" , which made it possible to 
perform the Brillouin zone (BZ) q-integration. On the 
other hand, the dynamical screening predicted by our 
LDA+DMFT calculations is detectable only over a short 
time scale. This essentially requires an integral over all 
frequencies, i.e., a cut-off flc 00 in 

2 _ 3 j"nc -Sbz X(q, w) b{^) dq duj 



= - lim / lmxioc{^) b{uj)duj. (6) 

Here b{x) = l/(e^^— 1) is the Bose-Einstein distribu- 
tion function, and the coefficient 3 comes (as mentioned 
also above) from the sum over the three spin components 
((52) = (5'2) + (5'2) + (5'2) in the PM phase, where we 
can consider the system to be magnetically isotropic). 

The most problematic step is the frequency integral 
since, e.g., in the case of CaFe2As2 the spin spectral func- 
tions have been measured only up flc ^ 60 — 80 meV. 
One can try to extrapolate the data to higher-frequencies: 
We consider here explicitly the values obtained^!^ for a 
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FIG. 4. (Color online) Cut-off dependence of the frequency 
integral of Eq. ((6| for the local (r = 0) and instantaneous 
[t — 0) magnetic moment calculated from the LDA+DMFT 
spin-susceptibility with (blue line: diagonal terms only; dark 
gray line: total) and without (light blue line) vertex correc- 
tions (i.e., just the bare bubble contribution). The most im- 
portant contribution to the frequency-integral comes from the 
energy region uj > 100 meV, i.e. from higher energies than 
the typical ones of the neutron experiments performed for the 
Fe-based superconductors. If considering a (rounded) typical 
cut-off of 100 meV (vertical red dashed line), it is evident that 
the estimate of Eq. ^ for the local and instantaneous mag- 
netic moment is much smaller (at least three times) than its 
"correct" asymptotic value. 



cut-off ric = lOOmeV, i.e., not too far away from the 
experimental window. 

With such a cut-off, the corresponding experimental 
estimate for mioc{t = 0) in CaFe2As2 obtained via Eq. 
ini (with an extra factor i for obtaining the squared mo- 
ment per Fe atom^) was 'mioc,t=o ^ 0.5/iB (to be com- 
pared with 0.95^B for the ordered moment in the same 
compound^). Notice that similar value for ?Ti/oc,i=Oj 
and namely OAfiB — O.S/is have been also estimated for 
BaFe2As2^i^ (for the optimally doped compound or in 
the ordered case, respectively). 

At first glance, hence, one would conclude that the 
size of the local and instantaneous moment in CaFe2As2 
is significantly smaller than that predicted by the 
LDA-I-DMFT calculations (note that for the 122 fam- 
ily, even larger value for mioc,t=o than for LaFeAsO are 
foundi^). However, one should keep in mind that the 
LDA-I-DMFT r = value corresponds to ilc = oo in 
Eq. dH). 

In order to understand to what extent the observed dis- 
crepancy might originate from such a cut-off "mismatch" 
between theory {Qc = +oo) and experiment (after a data 
extrapolationiii^ up to f^c = lOOmeV), we can now an- 
alyze further the frequency dependence of the local spin 
susceptibility, which was discussed in the previous sec- 
tion. From the bare bubble contribution to the spin sus- 
ceptibility (Fig. [31 left panel), it is clear that the major 



contribution to the integrated spectral weight comes from 
frequencies larger than leV. Including vertex corrections 
(Fig- El right panel) yields a richer spectral structure with 
a prominent peak at about u! 0.2eV. Our theoretical 
prediction of such a peak is not only consistent with the 
increase of q-integrated experimental neutron spectral 
fmiction for CaFe2As2 and Ba(Fei_2;Co2;)2As2 reported 
in Rcf. |48| ). but has also been very recently confirmed 
by novel neutron measurements^ in an extended energy 
range. 

In any case, the spectral weight lying below the ex- 
perimental cutoff represents just a small fraction of its 
total value, even the main features of the LDA-I-DMFT 
spectral function are located beyond the experimental 
cut-off. Hence we have plotted in Fig. |4] the value of the 
frequency integral over the LDA+DMFT data up to a 
cutoff flc- the experimental cut-off ftc ~ lOOmeV yields 
only 10% to 15% of the instantaneous squared local mo- 
ment Eq. (|6]). If assuming that also for CaFe2As2 only 
such a fraction of the moment has been integrated up to 
flc = lOOmeV, we get an estimate of mioc,t=o ^ l.S/is 
as the true magnetic moment. This represents of course 
only a rough estimate, significantly larger than the neu- 
tron spectroscopic values at about ftc = lOOmeV but still 
smaller than the LDA-^DMFT/XAS predictions. The 
rcmaing discrepancy might be partly originated by the 
lack of SU(2) symmetry of the interaction terms we as- 
sumed in Eq. ([l]), or, more hkely, by the effects of spatial 
correlations neglected by DMFT, which -at present- can- 
not be easily included via the cluster— or diagrammatic 
extension a^^'^^ of DMFT in complicate systems like the 
Fe-based superconductors. 



V. SUMMARY AND CONCLUSION 

In this paper, we studied by means of LDA+DMFT the 
local magnetic properties of iron-based superconductors 
focusing, in particular, on the 1111 systems (LaFeAsO) 
for which a small ordered magnetic moment (0.3 — 0.6/zb) 
is experimentally observed. We have constructed by 
the Wannicr-interpolation technique^i^ an effective low- 
energy model, which contains the dxy, d^^-yS, and dy^/zx 
orbital degrees of freedom, and supplemnted it with lo- 
cal interactions U = 1.8 eV and J = 0.45 eV. These in- 
teraction values arc motivated from constrained random 
phase approximation (cRPA) calculations and reproduce 
analogous results for the spectral function of the dpp 
modeli^i^: The system presents an intermediate value 
of the quasiparticle renormalization factor ^ 2 induced 
by electronic correlations, in agreement to the results of 
the ARPES experiments. 

In order to clarify the physics of local magnetic mo- 
ments and their dynamical screening as well as to make 
contact with the experiments, we have analyzed the lo- 
cal magnetic spin susceptibility within the framework of 
LDA+DMFT and looked into the detail of its tempera- 
ture and frequency dependence. As for the former, we 
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found the formation of considerably strong local (r = 0) 
and instantaneous (t = 0) magnetic moment, which is 
mainly driven by inter-orbital Hund's coupling J, and 
whose magnitude decays rapidly in imaginary time. More 
specifically, the characteristic features of such as "Hund's 
metal" physics are: (i) The instantaneously (short-time) 
local magnetic moment which corresponds to the square- 
root of the local spin correlation function xioc{t = 0) 
shows a saturation with decreasing temperature, (ii) 
The dynamically screened value at long times (t = /3/2), 
on the other hand, decreases upon lowering the tempera- 
ture, as a consequence of the still relatively good metallic 
properties of the system, [iii) The local magnetic sus- 
ceptibility at w = 0, which is the average (integration) 
over all r's, shows a (strongly renormalized) Pauli-like 
behavior in agreement with Rcf. [36l . For larger J values, 
a crossover to a Curie-like susceptibility is found^i. 

A clear hallmark of the physics of the dynamical 
screening of the large local instantanous magnetic mo- 
ments has been individuated in the magnetic spectral 
function, namely in the q-intcgrated spin susceptibilty 
Imx(w). In fact, when the vertex correction are prop- 
erly included, the total spin susceptibility calculated in 
LDA-I-DMFT diplays a characteristic, very pronounced 
peak around w ~ 0.2eV. 



When comparing our results to experiments, we note 
that the existence of a large (spin) magnetic moment 
on the Fe-atoms is supported by x-ray absorption spec- 
troscopy data for LaFeAsO^. The comparison with the 
neutrons spectroscopy data is more difficult: According 
to the LDA-^DMFT calculations about 90% of the mag- 
netic spectral weight lies above the typical experimental 
frequency window ^ 60—80 meV. This is compatible with 
the observed experimental increase of the q-integrated 
Imxioc(w) throughout the measured frequency range, as 
well as with the peak structure of the local spin suscep- 
tibility in one of the most recent neutron experiment^, 
and yields a crude estimate for the instantaneous mag- 
netic moment of at least ~ 1.5/is from the neutron data. 
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